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Abstract 

Following Sham and Rice [L. J. Sham, T. M. Rice, Phys. Rev. 144 (1966) 
708] the correlated motion of particle-hole pairs is studied, starting from the 
general two-particle Greens function. In this way we derive a matrix equa- 
tion for eigenvalues and wave functions, respectively, of the general type of 
collective excitation of a N-particle system. The interplay between excitons 
and plasmons is fully described by this new set of equations. As a by-product 
we obtain - at least a-posteriori - a justification for the use of the TDLDA for 

simple-metal clusters. 
36.40.+d, 31.50.+W, 33.20.Kf 
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I. INTRODUCTION 



Though there is no exact theorem sanctioning the use of time-dependent density function- 
als, the TDLDA (Time-Dependent Local Density Approximation) as intuitively introduced 
by Zangwill and Soven in 1980 Q turned out to be working extremely well in the descrip- 
tion of many-body effects in atoms PJ, molecules 0, clusters ||, and solids M. A common 
feature of all these examples is the charge density character of these collective states. But 
besides this type of correlated particle-hole motion there is another one well known from 
the optical properties of semiconductors, ionic solids, and rare gases as well, namely, the 
bound states of particles and holes, which generally occur if the attractive part of the general 
particle-hole irreducible interaction || dominates. In metals this does not happen, because 
the screened particle-hole interaction does not support bound states for usual values of the 
Thomas- Fermi screening vector ||. In small metal particles the situation is more subtle 
because all states are size-quantized 0]; this implies that the effective particle- hole inter- 
action is anywhere between a semiconductor (e 2 /er, with e the static dielectric constant) 
and a metal (e 2 exp(kr)/r, with k the Thomas- Fermi screening vector). For this reason, a 
more general theory is mandatory for the study of the optical absorption by collective states 
in metal clusters. Recalling the general discussion on the interplay between excitons and 
plasmons in the past ||[5| , the study of the general two-particle Greens function seems to 
constitute a promising tool to answer the question of whether or not the oversimplifying 
TDLDA has to be replaced by a more complete theory , in order to obtain a satisfactory 
insight in the nature of collective states in metal clusters. Starting from the method of Sham 
and Rice [|K]], a new effective equation is derived in section two, for the calculation of the op- 
tical spectra in metal clusters. Using the well-established jellium model J7|, these equations 
are solved for a number of Na-jellium clusters in section three. Finally, the conclusions and 
the comparison with alternative models of calculating the optical response in metal clustes 
are presented in section four. 
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II. THEORY 



We are interested in investigating the properties of excitons and plasmons in general 
many-electron systems. Therefore, our starting point is the full two-particle Green's func- 
tion, suitable to describe excitations of such a system involving the correlated motion of 
two quasi-particles. It will be shown that a unified description of excitonic and plasmonic 
excitations leads to a matrix equation which is formally similar to the well-known RPAE 
equations in atomic physics, with the essential difference of a better description of exchange, 
which becomes now properly screened. This formal similarity is, of course, not acciden- 
tal, and reflects the many possible ways one can derive the equations of motion for these 
excitations 01. 



Following Sham and Rice |IIJ we start by defining, field-theoretically, the exci- 
ton/plasmon wave-function, 

f x (fl,r 2 ,u) = < N,0\i>\f 2 ),1>(r-i)\N,\>= (1) 
= lim — f duxe^Go-xifl, Ui, r 2 , u - u x ) 

in terms of the appropriate function Gq ; \, which occurs in the Fourier transform of the 



general two-particle Green's function and is defined as follows |L0 



< N, 0Wfi)\N + l,a><JV + l, a\f(j$\N, A > 
G . x (ri, wi; r 2 , uj 2 ) = ^ — — - Q — (2) 

a LOi + 10 - h, N - & N+ i 

v < N, 0\^(r 2 ) \N - 1, a >< N - 1, a\ip(f{) \N, A > 

As usual, \N, > and \N, A > are the exact ground-state and plasmon/exciton excited state 
of the many electron system, while the tp are the usual field-operators. A represents a possible 
set of conserved quantum numbers. Gq;A satisfies the following homogeneous Bethe-Salpeter 



equation 10 



G 0;A (1, 2) = G(l, l')G(2', 2)7(1', 2'; 3', 4')G 0; a(4', 1') (3) 

where 1 = (f{,uji), etc. and integration over repeated indexes is assumed. I is the irreducible 
particle-hole interaction, which contains the essential ingredients of the physical processses 



we want to describe. For /, we consider the ladder approximation and neglect retardation 
effects, including the basic interactions described by the diagrams in fig. 1. Upward pointing 
arrows describe electron states, downward pointing arrows describe hole states, the dashed 
line represents the bare particle-hole interaction and the double dashed line represents the 
screened particle-hole interaction. As is well-known P,|TU|], since I is an irreducible interac- 
tion, only diagram b) must be screened (and vertex corrected). 
To proceed, we now explicitly write, for J, 

I(f{, r 2 , 1%, fl) = V b are(ri - f 2 )8{f 2 - r 3 )S(fl - r 4 ) - W(f{, f 2 )5{f{ - fz)8{f 2 - n), (4) 



where Vf, are is the bare interaction and W is the screened interaction ||11|| . Considering the 
usual expansion for f\ 



fx = X ph (j) p (j) h + Y ph (f) h *(j)p* 



(5) 



where X ph and Y ph are the forwards-going and backwards-going electron-hole amplitudes, 
respectively, whereas are quasi-particle wave-functions, and inserting this expansion, to- 
gether with the expression for I in eq.(2), and its solution in eq.(l) we obtain, after a final 
frequency integration, the following matrix-equation, 



(6) 
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where the matrices A and B are defined as 



A ph ,p'h' = (e p - e h )5 P:P >5 h> h>+ < ph'\I\p'h > 
B ph ,p>h> = < pp'\I\h'h> . 



(7) 
(8) 



In eqs.(7),(8), e p (e/i), are quasi-particle (quasi-hole) energies, whereas the matrix-elements 
between quasi-particle states read 

< ph\I\p'ti >= / dfl / dr 2 / dr 3 / dfl(f) p *{f{)())\{f 2 )I{f{, r 2 , r 3 , fl)(j) p/ {f^(j) hl {fl) (9) 
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Because these equations are formally identical to the RPAE equations, except for the residual 
particle-hole interaction, which has now a screened exchange term, we call the present results 
the RPA-SE equations. To our knowledge, this is the first time these equations have been 
explicitly obtained. Furthermore, we would like to point out that these equations are useful 
for the case of static screening, in which retardation effects are neglected [12]. However, the 



modifications to be introduced in the present formalism, in order to incorporate dynamical 
screening effects are formally straightforward. We defer the study of dynamical screening 
effects to a future publication. 

In order to implement the RPA-SE, we start from a independent particle picture. Tradi- 
tionally, the Hartree-Fock independent electron picture is the starting point. However, and 
since we are interested in studying the optical response of small clusters of simple metals, 
for which the Hartree-Fock approximation is known to provide a rather poor starting point, 
we start from a Density Functional Theory (DFT) scheme in the Local Density Approxima- 
tion (LDA). Therefore, we compute the one-electron wave-functions and energies in LDA to 
DFT, and correspondingly, the one-particle Green's functions. For the screened interaction 



W, we compute it in Hedin's GW approximation [r3|H|, following the prescription pro- 



posed by Hybertsen and Louie |15|. In this framework, instead of working out the complete 



self-consistent solution of the GW equations, one starts from the LDA approximation and 
carries out one iteration of the GW self-consistent scheme. This corresponds to write, for 
W, 



e 2 f e 2 e 2 
W(f[,f* 2 ) = 7-=, zrv + / df\df 2 — :7TX(r1,rl)— — , (10) 

\T\—T2\ J n — r 3 P4 — T2 



where, the density-density correlation function x is computed at the level of TDLDA. In 
general \ is a complex and energy-dependent function. In the static approximation, we take 
E — and \ becomes a real function. 

Within TDLDA, or equivalently, the matrix RPA-LDA (see below), the general form for 
the interaction / is replaced by the so-called residual interaction P,|rS|. It reads, according 
to the previous definitions, 



K(f{,fi,f^,fl) = 8(f 3 - r 2 )8(fl - f{)[V bare (f[, r 2 ) + Z1 ^- 8(fi - f 2 ) ] (11) 



o"( r ~3 - r!)5(rl - r~l)[ 



n - r 2 



+ —, — -r 2 ) . 
an 



III. RESULTS AND DISCUSSION 

In this section we apply the new theory to a number of magic jellium clusters of Na, and 
compare them to other models, mainly to the old fashioned TDLDA 0. Fig.2 gives the 
photoabsorption cross-section a per delocalized electron for various models of calculation. 
As a reference the dashed line gives for the results of the matrix-RPA-LDA. This curve 
is (of course) completely equivalent to the TDLDA. The only difference is that the latter 
calculates the dynamical polarizabilty and from this quantity the optical absorption, whereas 
the former calculates eigenvalues and eigenvectors, respectively, and from these quantities 



the optical absorption as described in detail by Yannouleas and Broglia [16|. That means, 
each spectral line is Lorentzian broadened (here with the parameter 77 chosen to be 0.1 eV). 

Next we solve the new equations, that means the irreducible particle-hole interaction is 
considered as given by eq.(4), meaning that the attractive part of the interaction is properly 
screened, and the repulsive part is calculated as usual. The result is given as a full line 
in fig.2. The difference to the matrix-RPA-LDA is not very spectular (in the specific case 
of Na%), because the exchange part of the interaction is less attractive. A little blue shift 
results for the main line of the optical response. 

Besides this effect there is a little redistribution of oscillator strength on the blue side- 
bands of the response on the high energy side of the main line. This can be made more 
visible by an alternative representation of the results in fig. 2a, which shows the line-spectrum 
(without Lorentzian broadening). That this is a general result can be seen in fig. 3, which 
shows our results for Na 2 o- In addition, the free-particle distribution of oscillator strength is 
given by crosses: upon inclusion of the electron-electron-interaction in the optical response 
the two strong lines which have been termed in the past as "fragmented plasmon " do not 
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change very much if one includes screening in the exchange part of the general interaction. 
As it was the case for Na 8 the main change occurs in the high-energy wing of the spectrum. 

As one can see in the next figures 4 and 5 for Na^s and Na$2 this is a general result: 
the proper screening of the exchange part is necessary on general reasons, but the detailed 
numerical consequences are rather negligible. In a sense this is a gratifying result, because it 
shows that still the simple TDLDA is the most efficient linear response-model for a variety 
of cases (atoms, clusters, surfaces and bulk). 



IV. CONCLUSIONS 

The main qualitative features of the optical response of metal clusters are given by the 
overwhelming importance of the formation of the colllective oscillation of all delocalized elec- 
trons: this phenomenon already present in the classical description of the optical properties 
H persists to exist in a trully microscopic description || using the rather simple model 
of the TDLDA; therefore it seems not too surprising that a true many-body solution of 
the problem, based on Feynman diagrams - in the way proposed by Sham and Rice |K| - 
does not change very much our picture of the optical response. But nevertheless a study 
along this line has to be done, because there is still an ongoing controversy about the right 



description of the effects of exchange in this context 

For this reason the marginal difference between TDLDA and the BSE description is 
highly welcome, because the latter is rather accurate but numerically complex, whereas the 
former is numerically easy to do but lacks a thorough justification. In this respect, we have 
a similar situation as we have with the so-called self- interaction correction SIC fi~8f , which 
has no sanctioning by any exact theorem, but turns out to be numerically equivalent to a 
quasiparticle calculation at the level of the GW-approximation [|11]] (which is considered to 
be on a firmer theoretical basis!). The general ph-interaction seems to be more important 
in clusters of other elements (e.g. rare gases, ionic species, or semiconductors), which will 
be studied next. 
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FIGURES 



Figure can be obtained from: ekardt@fhi-berlin.mpg.de 

FIG. 1. Graphical representation of I as used in this work. It originates from all parts of 
the irreducible particle-hole interaction of first order: bare Coulomb + bare exchange after 
the exchange part has been properly screened. According to the general rules for Feynman 
graphs part a) must NOT be screened, because this would result in a reducible diagram 
(compare ref. |10| ). Vertex corrections are not considered in this work. 
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FIGURE 2 



~i 1 r 



Na 



3 - 



8 

RPA-LDA 
" RPA-SE 



CM 
°< 



0)2 





FIG. 2. Photoabsorption cross section per delocalized electron for jellium Na$ in two 
different approximations: with a dashed line, the matrix- RPA-LDA fllfj, which is equivalent 
to TDLDA [[3]; with a continuous line, the RPA-SE as developed in this paper. The blue 
shift of the main line results from the proper screening of the attractive part of the general 
irreducible particle-hole interaction ( 110,11]). For details, see main text. 
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FIGURE 2a 
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FIG. 2a. Frequencies and oscillator strength of various particle-hole transtions in jellium 
Nag in three different approximations: crosses: the so-called free response (meaning, that 
effects of the induced charge are neglected), the rest is the same as in fig.l, but without 
Lorentzian broadening. We chose this way of representation, in order to make more visible 
the redistribution of oscillator strength in the high-frequency part of the spectrum. 
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FIGURE 3 



0.6 



0.4 



0.2 



- 


1 


1 1 1 

Na 20 ; 






[ . 1 FREE 






V//A RPA-LDA 






RPA-SE 











2 4 6 

E , eV 



FIG. 3. The same as fig. 2a, but for Na 2 o- 
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FIGURE 4 
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FIG. 4. The same as fig. 2a, but for Na 58 . This result is interesting, because it shows 
clearly a fragmentation of the line shape in the collective region around 3 eV. This has not 



been found in the preceding study of C. Yannouleas, et al. [19 
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FIGURE 5 
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FIG. 5. The same as fig. 2a, but for Na 92 . Microscopically, there is a fragmentation in 
the collective region. This, however, should be hard to resolve experimentally (which agrees 
with the interpretation of ref . JTIJ ) . 
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